function res = ExpandedChebyshev(ngp)
    res = zeros(ngp, 1);
    for k = 1 : ngp
        res(k) = cos((2*(ngp - k + 1.0) - 1.0)*pi/(2.0*ngp))/cos(pi/(2.0*ngp));
    end
end